Conductance of a disordered graphene superlattice 
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ON ' We study the conductance of disordered graphene superlattices with short-range structural cor- 

relations. The system consists of electron- and hole-doped graphenes of various thicknesses, which 
fluctuate randomly around their mean value. The effect of the randomness on the probability of 
transmission through the system of various sizes is studied. We show that in a disordered su- 
perlattice the quasiparticle that approaches the barrier interface almost perpendicularly transmits 
through the system. The conductivity of the finite-size system is computed and shown that the 
conductance vanishes when the sample size becomes very large, whereas for some specific structures 
the conductance tends to a nonzero value in the thermodynamics limit. 
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PACS numbers: 68.65.Cd, 73.22.-f, 73.63.-b, 73.40.Lq 



^3 ■ I. INTRODUCTION 

Graphene, a single atomic layer of graphite, has been successfully produced in experiment 1 -, which has resulted in 
intensive investigations on graphene-based structures, due to the fundamental physics interests that is involved and 
. the promising applications 2 -. There are significant current efforts devoted to growing graphene epitaxially 3 by thermal 
decomposition of silicon carbide (SiC), or by vapor deposition of hydrocarbons on catalytic metallic surfaces, which 
could later be etched away, leaving graphene on an insulating substrate. The low energy quasiparticle excitations 
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in graphene are linearly dispersing, and are described by Dirac cones at the edges of the first Brillouin zone. The 
linear energy- momentum dispersion has been confirmed by recent observations 4 . The slope of the linear relation 
i_i corresponds the Fermi velocity of chiral Dirac electrons in graphene, which plays an essential role in the Landau- 
\ Fermi liquid theory^ and has a direct connection to the experimental measurement. 

There are some unusual features of graphene, such as the effects of electron-electron interactions on the ground-state 
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I properties^, anomalous tunneling effect described by the Klein tunneling, the tunneling through a p-n junction^ that 
follows from chiral band states, and the energy-momentum linear dispersion relation. The Klein tunneling predicts 
that the chiral massless carrier can pass through a high electrostatic potential barrier with probability one, regardless 
of the height and width of the barrier at normal incidence, which is in contrast with the conventional nonrelativistic 



massive carrier tunneling where the transmission probability decays exponentially with the increasing of the barrier 
hight and would depend on the profile of the barrier iV ' 11 ! 12 

An exciting experimental development is the ability to apply an electric field effect or submicron gate voltage, 
in order to illustrate graphene p-n junctions* 1 - 3 - By applying an external gate voltage, the system can be switched 



from the n-type to the p-type carriers, thereby controlling the electronic properties that give rise to graphene-based 
nanodevices. Recently, strong evidence for Klein tunneling across potential steps which is steep enough in graphene 
has been experimentally observed^ 

Clean graphene junctions were predicted to display a number of fascinating physical phenomena, even in the absence 
of electron-electron interactions^ Interestingly, the Veselago lensing of electric current by a single p-n junction in 
clean graphene 1 -^ and the Andreev reflection, and the electron to hole conversion at the interface at normal incidence 1 ^, 
have all been predicted theoretically. Such phenomena are predicted to change both quantitatively and qualitatively 
when disorder is included in the model. For instance, inhomogeneous graphene p-n junction systems were studied 
using the Thomas-Fermi approximation, including disorder effects, by Fogler and collaborators^. They showed that 
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junction resistance is dominated by either ballistic or diffusive contributions, depending on the density of charged 
impurity and gradient of the carrier density. 

In the semiconductor context there are basically a large number of works on the tunneling, which have resulted in the 
"obvious" declaration that the electronic properties of semiconductor supcrlattice are different from those calculated 
in a single-barrier junction. Moreover, the electronic properties of semiconductor superlattices in the presence of 
disorder have been studied by several group o 19 ' 20 i 21 ' 22 i 23 . Importantly, all the electronic states are localized in the 
thermodynamic limit for a semiconductor superlattice in the presence of white-noise disorder. 21 

Graphene superlattices, on the other hand, may be fabricated by adsorbing adatoms on graphene surface through 
similar techniques, by positioning and aligning impurities with scanning tunneling microscopy 2 ^ or by applying a 
local top gate voltage to graphene^ Recently, a periodic pattern in the scanning tunneling microscope image has 
been demonstrated on a graphene on top of a metallic substrate^ The transition of hitting massless particles in 
graphene-based superlattice structure (GSLs) was first studied by Bai and Zhangi^l They showed that the conduc- 
tivity of the GSLs depends on the superlattice structural parameters. Furthermore, the superlattice structure of 
graphene nanoribbons has been recently studied by using first-principles density functional theory calculations. 28 
These calculations showed that the magnetic ground state of the constituent ribbons, the symmetry of the junction, 
and their functionalization by adatoms represent structural parameters to the electronic and magnetic properties of 
such structures. Recently, novel physical properties of GSLs with one-dimensional (ID) Kronig-Penney type and 
2D muffin-tin type potentials were also studied^ The results showed that a periodic potential applied by suitably 
patterned modifications leads to further charge carrier behavior. The propagation of charge carriers through such 
a superlattice is highly anisotropic, and in extreme cases results in group velocities that are reduced to zero in one 
direction but are unchanged in the other direction. Moreover, they showed that the density and type of carrier are 
extremely sensitive to the applied potential. 

It would, therefore, be worthwhile to investigate how the conductance of graphene superlattice junctions are af- 
fected by structural white noise, and compare the conductances with those calculated for disordered semiconductor 
superlattice. Due to the conservation of pseudospins in graphene, backscattering process is suppressed at normal 
incidence, which makes the disordered regions transparent. 1 

The purpose of this paper is to study the electronic behavior of graphene superlattices p-n junctions by using 
the transfer-matrix method. The system that we study consists of a sequence of electron-doped graphene as wells, 
and hole-doped graphene as barriers. We study the effect of the disorder imposed on the size of the barriers in the 
transmission probability, T, through the system as a function of the system size (number of the barriers), together 
with various incident angles. The dc conductance of the finite-size system takes on a nonzero value of the transmission 
in some special configurations. Using the finite-size scaling of transmission, we show that the conductance, in the 
thermodynamic limit, tends to a finite constant for spacial cases. 

The rest of this paper is organized as follows. In Sec. II we introduce the models and derive the related transfer 
matrix. We also explain how we calculate the transmission probability and the dc conductivity. Section III contains 
our numerical calculations. We conclude in Sec. IV with a brief summary. 

II. MODEL AND THEORY 

Consider a system of superlattice p-n junctions in the independent carriers model at zero temperature, and in the 
absence of carrier-phonon and spin-orbit interactions. The low-energy massless Dirac-band Hamiltonian of graphene 
in the continuum model can be written a o 30 ' 31 

Ho = hvr (<r 1 ki + cr 2 k 2 ) 
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where r = ±1 for the inequivalent K and K 1 valleys at which it and it* bands touch, fcj is an envelope function 
momentum operator, v is the Fermi velocity, and Oi is a Pauli matrix that acts on the sublattice pseudospin degree 
of freedom. The total Hamiltonian of a massless carrier in a special geometry is written as, H — Hq + V(x) where 
V(x) is the graphene-based superlattice potential which is modeled as described below. 



A. Superlattice Model 



We consider superlattice p-n junctions in a graphene-based structure. The system consists of two kinds of graphene 
with different potentials, the first being an electron-doped graphene with thickness dyy, while the second is a hole- 
doped part with thickness ds, standing alternately. The potential for the electron- and hole-doped graphene are Vq 
and zero, respectively. The energy of the incident particle is Eq = 2irhv/A with the wavelength A across the barriers, 
in such a way that the Fermi level lies in the conduction band outside the barrier and the valence band inside it, i.e., 
(0 < Eq < Vq), as shown in Fig. 1. The growth direction is taken to be the x axis which is designed as the superlattice 
axis. In order to neglect the strip edges, we assume that the width of the graphene strip is much larger than ds- We 
set disorder situations in which the value of ds fluctuates around its mean value, given by < ds > = b. In the model 
the fluctuations are given by, ds\i = b(\ + 8 u), where {e^} is a set of uncorrelated random variables or white noise 
with the box distribution, — 1 < e< < 1, and i is the site index. Here, the 5 is the disorder strength. 
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FIG. 1: Model of graphene superlattice p-n junctions. 



We consider graphene-based superlattice potential in a simple model as 

dn\i 



V(x) = \ Vq if \ x ~ X2l \ < 
otherwise , 



(1) 



where X2i is the position of barriers' center. The model is similar to the potential of semiconductor superlattices that 
has been used by other groups.— 



B. DC conductivity 



Let us now consider the case in which the incident massless electron in the GSLs propagates at angle along the 
x axis (see Fig. 1) and, therefore, the Dirac spinor components, "01 an d "02, which are the solutions to the Dirac 
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Hamiltonian, can be expressed^ as: 



M*>V) = (aie iKiX + b ie - iK * x )e ik » y 

M*,V) = Si{aie iKiX+i ^ -b i e- iKiX - i * i )e ik vV , (2) 



where 

E 



Si = sgn(E - V(x)), k y = — sin(^) , (3) 

nv 

and 

{k x = Eq cos(<j))/hv for well 

q x = yj (Eq — Vo) 2 /H 2 v 2 — k 2 for barrier 

In order to calculate the transmission coefficients, we use the transfer-matrix method^. To this end, we apply the 
continuity of the wave function at the boundaries, and construct the transfer matrices as follows 

r ) ~ 2^4, \ e» + e ie e» - e' ie ) F{ZN) X (5) 

e lk ^( e ~' e - e ^)/[2e <g ^« cos6»] \ 
e lk ^ (e te + e^)/[2e- i9 »'» cos6»] ) 2N 

where r and t2N are the reflection and transmission coefficients of the system that consists of N barriers, and p(2N) 
is the transfer matrix given by 

2N 

P(2N)= Y[Pi,i-i (6) 

«=3 

Mu M12 



Pi 



2,2—1 



M 2i M 22 

where also P^j-i is a transfer matrix from site i to % — 1 and My are given by 

Mu = e iKil »-^e- iK ^ l ^[e- iv ^ - e i * 3 *]/2cos(^ (i _i ) ) (7) 

M12 = e^^'^-De-^-D't-Dle-^t- 1 ) + e -^]/2 cos(^ (i _i)) 

M 2 i = e^^-V^- 1 ^- 1 '!^- 1 ) +e iv *]/2cos(^ (i _i 3 ) 

M 22 = e -«^i'(*-i) e «(*-i) '(*-!) [ e ^«-i) -e-^J^cos^^i)) 

where = Xj*"* 1 *^ d B |j + »nt[(i - l)/2] dw is the length of system at i-th boundary and moreover, 

{</> for well ,_, 

= tan _1 (fcj,/q :E ) for barrier * ' 

It is evident that T(Eq, 4>) = |*2iv | 2 , and that it can be calculated from Eq. ([5]) for a given N. When the transmission 
coefficients are calculated, the conductivity of system is computed by means of the Biittiker formula^ taking the 
integral of T(E , <j>) over the angle 



G = G T(E, 4>) cos(»d0 (9) 
where Go = e 2 mvw/H 2 with w being the width of the graphene strip along the y direction. 
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III. RESULTS AND DISCUSSION 

Let us first calculate the transmission probability and study the electronic properties of disordered GSLs as a 
function of the strength of disorder introduced in the system. We consider the width of barriers as a random variable, 
so that the length of system in the numerical calculations will be L = N(b + dyy). In all the numerical calculations, 
we assumed b =< ds >= 50 ran, while the wavelength of the incident particle is set by A = 50 nm, or, equivalently, 
the energy of the carrier, Eq = 83 meV. In all of the calculations we used, Vq = 200 meV, unless otherwise specified. 
The number of realization of the random configurations is about 500. 

Figure 2 shows the transmission probability, T, of the incident electrons hitting a GSLs, as a function of the angle <f> 
for several values of the disorder strength, S. The number of the barriers in the figure is, N — 100, with dw = 10 nm. 
It is clear that, the transmission decreases by increasing the disorder for all the angles, apart from the strictly normal 
incidence case, = 0. This is physically understandable due to the Klein tunneling process in graphene, where the 
backscattering process is suppressed. Moreover, the massless carriers with incident angle close to normal incidence 
can survive in the presence of disorder, while the width of the angles around the normal incidence decreases with 
increasing strength of the disorder as well. 

In order to understand the finite-size effect and the effect of the width of the wells, the transmissions probability 
of a massless particle through the system were calculated as a function of the incident angle. The results are shown 
in Fig. 3 for several sizes at S = 0.1. In Fig. 3(a) we set dw — 10 nm. The transmission decreases with increasing 
system size for all the angels, except again at cj> = 0. This behavior is in contrast with a clean GSLs result, where 
the number of the peaks increases by increasing the number of the barriers.— In Fig. 3(b) the width of the wells is 
set to, d w = 30 nm. Two sharp peaks in the transmission are obtained that disappear in Fig. 3(a). Furthermore, 
in the case with dw — 50 nm there is only one angle, </> ~ 60°, that the transmission survives in the presence of the 
disorder, whereas under other angles the transmissions are suppressed. Meanwhile, there is clearly a wide domain 
around <j) — for which the transmission survives, and is larger than the one shown in Fig. 3(b). The width of the 
domain decreases with increasing strength of the disorder. It is worthwhile to note qualitatively that for 5 = 0, when 
we have the condition that, (q x + k x ) x (b + dw) = 2mir ( m is an integer), the transmission has finite values at angles 
different from <fi = 0. This is due to the resonance process in a system with N barriers. 

We also studied how disorder, introduced in the GSLs, affects the conductivity of the system. Hence, we also 
calculated numerically the dc conductivity by using Eq. ([9]), with the white- noise structural disorder imposed on the 
system. Figure 4 shows the dc conductance of the GSLs as a function of Vq for various strengths of the disorder, S. As 
shown in the figure, the conductivity of the GSLs decreases by increasing the strength of the disorder. However, the 
conductivity approaches a finite value, i.e., the existence of a finite conductivity in finite-size disordered GSLs should 
be expected. In general, the resonance condition is given by a function that yields, f(q x , dw, ds) — rem. For instance, 
for the case N = 1, the condition yields, q x ds = rrm, as a result of which T(<fi) would be an oscillating function of q x . 
Note that q x is determined by Vq . Consequently, this leads to a finite dc conductivity which is an oscillating function 
of Vq. The observation of conductance oscillations in extremely narrow graphene hetrostructures has been observed 
experimentally. 34 

In Fig. 5 the conductivity of the system is plotted as a function of Vq, where the width of the wells is d w = 10 
nm with S = 0.1. The conductivity decreases with increasing the size of system. The inset in the figure shows the 
conductivity of a clean GSLs as a function of Vq for several system sizes. It indicates that the dc conductance of 
clean superlattice behaves uniquely for different sizes, but in a disordered GSLs it decreases by increasing the size of 
system, as shown in Fig. 5. At a constant strength of the disorder, changing dw may also change the conductivity, 
as depicted in Fig. 6. It demonstrates that the conductivity varies periodically with increasing dw- As a result, in 
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the disordered GSLs, the dc conductance of finite-size systems depends on the structural parameters, especially dw- 
To compute all results that have been presented so far, we considered a system of finite size. Next, we wish to 
calculate the finite-size scaling of G/Gq. For this purpose, we calculated the conductivity as a function of the system 
size. The results are summarized in Fig. 7. Importantly, the conductivity vanishes by a simple power law, except 
for the case for which, A = dw = 50 nm. In general, for A = mdw the conductivity approaches a finite value as N 
becomes large. 

In order to examine such results better, we also calculated the G/Gq for a case for which A = dw = 45 nm. We 
found that the conductance tends to a constant in the thermodynamic limit, the numerical data are fitted by using 

|=^ + ^ ^ 

where g^, 7, and r\ are constants, g^ is the asymptotic value of G/Gq in the thermodynamic limit, N — > 00. As a 
result, for the case dw = 5 nm we obtain g x = 0, 7 ~ 1.0, and 77 ~ 0.46; for dw — 10 nm we obtained, j M = 0, 
7 ~ 0.9, and rj ~ 0.42, and g x ~ 0.14, 7 ~ 0.6 and 77 ~ 0.2 for dw = 50 nm. In all the case the regression was 
with r 2 — 0.99, indicating very accurate fits. Note that g^ is zero for small dw, but tends to a nonzero constant for 
dw = 50 nm. 



IV. CONCLUSION 



We studied numerically the dc conductance of a discorded graphene supcrlattice p-n junctions for various values 
of the strength of structural disorder imposed on the material. It was shown that there exists a width around the 
normal incidence angle for which the transmission becomes finite in the presence of structural white-noise disorder. 
That is, the white-noise disorder gives rise the largest number of the peaks in the transmission, suppressed in the 
thermodynamic limit but quasiparticles which approach almost perpendicularly to the barriers transmit through the 
material. We also calculated the conductivity of a finite-size disordered system and showed that the conductivity 
decreases by increasing the system size but that there are cases for which the conductance approaches a nonzero 
value. This result is in contrast with the case of a clean (ordered) GSLs.— Furthermore, the results of the finite- 
size scaling computations predict a zero conductance for all the GSLs, except for some special dw values for which 
A = mdw, where m is an integer, in which case the conductance tends to a nonzero constant in the thermodynamic 
limit. 

Apparently, such a feature is independent of the value of b =< ds >■ Consequently, we predict a finite conductivity 
for a disordered GSLs when the wavelength of incident particle is equal to mdw- These results are in complete 
contrast with those calculated for disordered semiconductor superlattice which become insulator i 19 i 21 ' 22 Our finding 
for the dc conductance of the GSLs should be important to the design of electronic nano-devices based on graphene 
superlattices. It would probably worthwhile to extend the present work to the case in which a correlated noise is used. 
In this case one must replace the white-noise with a proper short- or long-range correlated noise. 
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FIG. 3: (Color online) Transmission probability T of the massless carriers through the system as a function of the incident 
angle, the system size and 8 = 0.1 for dw = 10 nm. (a) dw = 30 nm. (b) dw = 50 nm. 
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FIG. 4: (Color online) DC conductivity as a function of the barrier potential, Vo (in units of meV) for various strengths of the 
disorder and N = 100 and dw = 10 nm. 




FIG. 5: (Color online) DC conductivity as a function of the barrier potential Vo (in units of meV) and system size for S — 0.1 
and dw = 10 nm. Inset shows the same, but for clean (ordered) GSLs. 




FIG. 7: (Color online) Finite-size scaling of the dc conductivity as a function of size N for the disorder strength 5 — 0.1 and 
various values of dw and A (in the inset). 



